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The effects of stellar population synthesis on the distributions of the 
asteroseismic observables Umax and Az/ of red-clump stars 
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ABSTRACT 

The distributions of the frequencies of the maximum oscillation power (I'max) and the large 
frequency separation (Aj/) of red giant stars observed by CoRoT have a dominant peak. 
Miglio et al. (2009) identified that the stars are red-clump stars. Using stellar population syn- 
thesis method, we studied the effects of Reimers mass loss, binary interactions, star formation 
rate and the mixing-length parameter on the distributions of tlie I'max and Ai/ of red-clump 
stars. The Reimers mass loss can result in an increase in the i^max and Ai/ of old population 
which lost a considerable amount of mass. However, it leads to a small decrease in those 
of middle-age population which lost a little bit of mass. Furtliermore, a high mass-loss rate 
impedes the low-mass and low-metal stars evolving into core-helium burning (CHeB) stage. 
Both Reimers mass loss and star formation rate mainly affect the number of CHeB stars with 
i^max and Ai^, but hardly affect the peak locations of i'max and Ai^. Binary interactions also 
can lead to an increase or decrease in the Vmax and At/ of some stars. However, the fraction 
of CHeB stars undergoing binary interactions is very small in our simulations. Therefore, the 
peak locations are also not affected by binary interactions. The non-uniform distributions of 
i^max and Av are mainly caused by the most of red-clump stars having an approximate radius 
rather than mass. The radius of red-clump stars decreases with increasing the mixing-length 
parameter. The peak locations of lymax and At/ can, thus, be affected by the mixing-length 
parameter. 
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1 INTRODUCTION 



Solar-like oscillations have been con finned in many main -sequence 
and subgiant stars, such as q Cen A JBedding et alj2004[). a Cen B 
I Kieldsen. Bedding & Butleil2005l). Procyon A JEggenberger et al.l 



l2004h . n Bootis Carrier. E ggenberger & Bouchv' '20051). etc. 
Some giant stars, s uch as ^ Hya ( Frandsen et al. 2002), r? Ser 
JBarban et al.l l2004h . and e Oph JDe Ridder. Barban & Carrier! 
|2006|) . have also been found the behaviour of oscillating. Astero- 
seismology is a powerful to ol for determining the fundamental pa- 
rame ters of individual stars ( Eggenberger, Carrier & Bouchyl 20051 : 
lEggen b erger & Carrier 2006 ; Yang & Meng 2010) , and it has sig- 
nificantly advanced the theories of stellar structure and stellar 
evolution. 

The detection of solar-like oscillations in red giant stars 
opened up a new field that can be explored with asteroseis- 
mic techniques. For an ensemble of cluster stars all having the 
same age and metallicity, studying oscillations in the cluster stars 



should provide more constraints on the theory of stell ar evolution 

than j ust fitting parameters of single oscillating stars iStello et al.l 
I2007D . Thus a number of attempts have been made to detect 
solar- l ike o scillations 
19931: loiuiland 
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in giant stars of c lusters ( Gilliland et al 
2008; Edmon ds & Gilliland 1996: Stello et al. 



2007 ; iFrandsen. Bruntt GrundaK 



Stello etal J i201Q : iBedding et al 



20071 : IStello & Gilliland l2009l : 
The oscillations of gi- 
ant stars have been detected in a few clusters, such as open 
cluster M67 JSt ello et al. 2007), globular cl uster 47 Tucanae 
JEdmonds & Gilliland .1996.) and NGC 6397 (iStello & GilUland 
120091) . However, these observations did not obtain oscillation fre- 
quencies. Rece ntly , using the data observed by the first CoRoT 
(Baglin et al. 20061) 150-day long run in the direction of the galac- 
tic centre (LRcOl). iHekker et al.l J200 9) obtained the frequencies 
of maximum oscillation power (Vmax) and the large separations 
(Av) of about 800 solar-like oscillating red giants. They found 
that distributions of the Vmax and Av are non-uniform (see Fig. 
m, which provide an opportunity for detailed studie s of popu- 
lation of galactic-disk red giants JMiglio et al J 120091) . Addition- 
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Figure 1. The histogram of the Umax and At/ obtained by iHekker et alJ 
i2009h . 



ally, theoretical Umax and A;/ can be calculated using equations 
JKieldsen & Beddindl 19951) 
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(2) 



and 

Thus, we can obtain the theoretical distributions of Vmax and Ai^ 
by simulating the distributions of stellar mass, radius and effective 
temperature. 

Using t he ste llar population synthesis code TRILEGAL, 
iMiglio et alj ( 120091) identified that the oscillating giants observed 
by CoRoT are primarily red-clump stars, i.e. post-flash core-He- 
burning (CHeB) stars, and found that theoretical distributions of 
Vmax and ^v are good global agreement with observed ones. 
Equations [T] and |2] show that the theoretical distributions of Vmax 
and Ai/ are relative to the mass and radius distributions of popula- 
tion. Therefore, Miglio et al. ( 2009) suggested that mass-loss rate 
during the red-giant branch (RGB) and star formation rate may af- 
fect the distributions of the Vmax and At/. 

In this paper, we focused the effects of mass loss, binary 
interactions, the mixing-length parameter, etc. on the distribu- 
tions of Vmax and Ai/ of CHeB s tars. We used the Hurley 
rapid single and binary evo lution codes JHurlev. Pols. & Toull200d : 
iHurlev, Tout & Polsll2002l) to construct stellar models. In Hurley's 
codes, mass-loss efficiency is an adjustable parameter. They can, 
thus, be applied to studying the effect of mass loss on the distri- 
butions. 

The paper is organized as follows. We simply show our stel- 
lar population synthesis method in section 2. The results are rep- 
resented in section 3. In section 4, we discuss and summarize the 
results. 



2 STELLAR POPULATION SYNTHESIS 

To simulate the s tellar population observed Vmax and Ai^ 
JHekker et al J 120090 . we calculated single-star stellar population 
(SSP) and binary-star stellar population (BSP), respectively. Stellar 
samples are generated by the Monte Carlo simulation. The basic as- 
sumptions for the simulations are as follows, (i) Star formation rate 
(SFR) is assumed to be a constant ov er the past 15 Gyr . (ii) Th e 
age-metallicity relation is taken from Rocha-Pinto et al. ( 2000a). 
(iii) The lognormal initial mass function (IMF) of Chabrier (2001) 
is adopted. Additionally, for BSP, we generate the mass of the pri- 
mary, Ml , according to the IMF. The ratio (g) of the mass of the 
secondary to that of the primary is assumed to be an uniform dis- 
tribution within 0-1 for simplicity. Then the mass of the secondary 



star is given by qM\ . We assume that all stars are members of bi- 
nary systems and that the distribution of separations is constant in 
log a for wide binaries and falls off smoothly at close separation: 



ania) 



{a/aoY 



ao < a < ai, 



(3) 



where Q^cp ~ 0.070, a,, = WRq, ai = 5.75 x IO^Rq = O.lSpc 
and m ~ 1.2. This distribution implies that the number of wide 
binary system per logarithmic interval is equal, and that approxi- 
mately 50% of the stellar s ystems are binary systems with orbita l 
periods less than 100 yr (Han, Podsi adlowski. & Eggletonlll995[) . 
With these assumptions, we calculated the evolutions of 10^ stars 
with an initial mass between 0.8 and 5.8 Mq and mixing-length pa- 
rameter a = 2.0. The theoretical Vmax and Av of single and binary 
stars are calculated using equations (O and ([2J. 



3 CALCULATION RESULTS 

3.1 Mass-loss effect 

Fig. [2] shows the distributions of the mass and radi us of the CHeB 
stars of simulated SSP with different Reimers JReimersI 119751) 
mass-loss efficiency (r?). The distributions of the initial mass of the 
progenitors of the CHeB stars also are shown in Fig. |2] Reimers 
mass loss mainly affects the stars with initial mass less than 2 Mq. 
With the increase in the ry, the mass lost by a star during the RGB 
increases, which leads to many low-mass stars can not evolve into 
the CHeB stage, and the lower limit of the initial mass of stars 
being able to evolve into the CHeB stage increases too. For exam- 
ple, when the value of 77 increases from 1.0 to 2.0, the lower limit 
increases from about 0.95 Mq to 1.25 AIq. Thus the fraction of 
low-mass CHeB stars and the sample size of the simulated CHeB 
stars decreases with increasing ry. For 77 = 2.0, the mass distribution 
is almost uniform. Moreover, the distribution of logarithmic radius 
has a fixed peak location between 1.05 and 1.15, which is almost 
not affected by the mass loss. 

The distributions of the Vmax and Av of the simulated CHeB 
stars are represented in Fig. [3] Both distributions have a dominant 
peak. But the peak locations are slightly lower than the observed 
ones. With the increase in 77, the number of the CHeB stars with 
Vmax in the range of about 10~30 fiUz and Au in the range of 
about 1~4 fiHz decreases. However, the dominant peak for Au al- 
ways locates between 3 /xHz and 3.5 ^Hz, which implies that the 
effect of the mass loss on the Av is not enough to change the peak 
location. For r; ^ 1, the dominant peak of Vmax locates between 
25 fiUz and 30 /iHz. However, when the r; increases from 1.0 to 
2.0, the peak location of Vmax moves to between 20 ^Hz and 25 
^Hz. 

The changes in the distributions of Vmax and Av shown in 
Fig. |3] should not be due to the typical variations in Monte Carlo 
simulations. In order to show the effect of sample size on the distri- 
butions of Vmax and Av, we calculated the evolutions of 10** stars 
which produced 3973 CHeB stars and the evolutions of 2x10^ 
stars which given 793 CHeB stars. Both samples have the same 
evolutionary parameters. Figure |4] shows that the distributions of 
Vmax and Av of the two samples are very similar. A Kolmogorov- 
Smirnov test shows that the discrepancies between the two sample 
distributions are not significant: the Kolmogorov-Smimov statistic 
D ~ 0.030 for Vmax and 0.033 for Ai/, and the significance level 
of the D is 0.61 for Vmax and 0.45 for Au. 

In Fig.|5] we plotted the distributions of the Vmax and Au of 
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Figure 2. The histograms of the mass and radius of the CHeB stars of sim- 
ulated SSP with different mass-loss coefficient rj. In the left panels, empty 
bars illustrate the mass distribution of the progenitors of the CHeB stars. 



the simulated CHeB stars with different 77 as a function of stellar 
age. For age > 2 Gyr, Fig.[5]shows that the value of the i^max and 
Av of CHeB stars is mainly located in the range of '--^ 10 - ~30 fiRz 
and ^1 - ~4 /iHz, respectively. Moreover, most stars are located 
around the upper boundary of the distributions. These stars are 
close to the zero-age horizontal branch (ZAHB) in the Hertzsprung- 
Russell (HR) diagram. While the stars around the lower boundary 
are approaching the asymptotic giant branch (AGB). The CHeB 
stars with higher Umax and Av are mainly from the young popula- 
tion with age less than 2 Gyr. Additionally, Fig.|5]clearly indicates 
that Reimers mass loss mainly affects the Umax and Az^ of old stars. 
A big mass-loss rate impedes the low-mass low-metal stars evolv- 
ing into CHeB stage. Consequently, for 77 = 0.5, 1.0 and 2.0, the 
fraction of the CHeB stars with age larger than 12, 10 and 5 Gyr 
is negligible, respectively. However, the mass loss scarcely affect 
the distributions of the Umax and Au of stars with age < 2 Gyr. 
When ri increases from 0.5 to 1.0, the upper boundary of the Umax 
and Au of the CHeB stars in the age range of about 7-10 Gyr is 
enhanced, but the effect of the increase in the rj on the distributions 
of the Umax and Au of the stars in the age range of around 2-7 
Gyr is not significant. When 77 increases to 2.0, the upper bound- 
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Figure 3. The histograms of the Umax and Au of the CHeB star's of simu- 
lated SSP with different mass-loss coefficient rj. 
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Figure 4. Histogram showing the effect of sample size on the distributions 
of the Umax and Au of simulated CHeB stars. Upper panels show the re- 
sults of 10'' model evolutions. While lower panels represent the results of 
2x10^ model evolutions. Both samples have the same other parameters. 
A K-S test shows that the discrepancies between two distributions are in- 
significant. 



4 W. Yang, X. Meng and Z. Li 






N O - 

X O - 



7; = 0.5 



10 15 




:l^ 



w 



5 10 15 

Stellar age (Gyr) 



5 10 15 

Stellar age (Gyr) 



Figure 5. The distributions of the Umax and Au of simulated CHeB stars 
with different »7 as a function of age. 



ary of the Vmax and Av of the stars with age approximately 5 Gyr 
Increases too, however, the distributions of the Umax and Au of 
the stars in the age range of about 2-4.5 Gyr move down slightly, 
which leads to the peak location in the Umax histogram moving 
to 20-25 ^Hz. 



3.2 BSP effect 

Fig.[6]shows the histograms of the Umax and Au of the CHeB stars 
of simulated BSP with r; = 0.5 and the distributions of the Umax and 
Av as a function of stellar age. The distributions of the Umax and 
Ai^ are very similar to those of the SSP with the same rj. This is 
because that most binary stars evolved into CHeB stage are wide 
binary stars. Binary interactions such as mass transfer and mass 
accretion do not take effect in these wide binary stars. However, 
the lower panels of Fig. |6]show that the value of Umax and Au 
of some stars with age > 2 Gyr is larger than 40 /^Hz and 4 /^Hz, 
respectively. This is due to the fact that binary interactions lead to 
an obvious increase or decrease in the mass of the stars when they 
evolve into the CHeB stage. Consequently, the Umax and Au of 
these stars are larger than those of the wide binary stars. However, 
the fraction of these stars is very small. Thus the effect of binary 
interactions on the distributions of Umax and Au is not significant. 



3.3 SFR effect 

The star formation rate of the Galaxy is not a constant. 
iRocha-Pinto et alj ( l2000bl) gave that the SFR of the Galaxy at 2- 
5 Gyr ago is about 2 times larger than an average SFR over the past 
15 Gyr and the SFR of the Galaxy at 7-9 Gyr ago is around 1.5 
times larger than the average SFR. Fig.|7]shows the distributions of 
Umax and Au of CHeB stars of the si mulated SSP with a constan t 
SFR and with the Galaxy SFR given bv lRocha-Pinto et al.1 i2000bh . 
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Figure 6. Upper panels show the histograms of the Umax and Au of the 
CHeB stars of simulated BSP. Here, rj = 0.5. Lower panels indicate the 
distributions of the Umax and Au as a function of stellar age. 
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Figure 7. The histograms of the Umax and At' of simulated CHeB stars for 
different star formation rate. Here, r] = 0.5. Upper panels show the results 
of SSP with constant SFR. W hile lower panels re present the results of SSP 
with the Galaxy SFR given bv iRocha-Pinto et al.i (2000b) . 



Although the number of simulated CHeB stars is changed, the dis- 
tributions are similar. In fact, for rj = 0.5, Fig.|5]shows that the value 
of the Umax and Au of the CHeB stars with age > 2 Gyr is mainly 
located in the range of 10-30 ^Hz and 1-4 fiKz, and gathers around 
25 /iHz and 3 fjHz, respectively. Thus increasing or decreasing the 
SFR of stars at a certain age > 2 Gyr can increase or decrease the 
number of CHeB stars but cannot affect the peak locations of the 
Umax and Au of CHeB stars. Although increasing the SFR of the 
young stars in the age range of 0-2 Gyr also can increase the num- 
ber of stars with relatively high Umax and Au, the peak locations 
cannot also be affected unless the SFR is enhanced many times. 



3.4 Mixing-length parameter effect 

The distribution of the radius and of the Aj^ of simulated CHeB 
stars have a dominant peak, but the distribution of the mass of the 
simulated CHeB stars with rj = 2.0 has not a dominant peak, which 
implies that the non-uniform distributions of Umax and Au may be 



The ejfects ofSPS on the distributions of the asteroseismic observables i^max <^nd Au 




Mass (M„) 





o 
o 
o 

o 
o 
in 



77 = 0,5 a = 2.4 



m^f^ 



50 100 

V (uHz) 



2 4 6 
t\v (/^Hz) 



10 



Figure 8. The histograms of the u^ax, Ai/, M, and R of the CHeB stars of 
simulated SSP with the mixing-length parameter a = 2.4. 



caused by the radius of many CHeB stars concentrating in a narrow 
interval. In fact, equations [T] and [2] indicate that Vmax and Au are 
more sensitive to the change in radius than that in mass. Further- 
more, a variation of the mixing-length parameter a mainly changes 
the stellar radius, but has alm ost no influence on the luminosity 
dKip penhahn & Weieert 1990). Used the Eggleton's stellar evolu- 
tion code CEggleton. 197 L . 1972. . . 1973.) . our calculations show that 
the increase in a leads to a decrease in the radius of stars evolved 
into horizontal branch, and thus leads to an increase in the i^max 
and Ai/of the stars. In order to study the effect of the mixing-length 
parameter a on the distributions of Umax and Au, we modified the 
Hurley's single evolution code to calculate the stellar evolution with 
a = 2.4. Fig.[8]shows the histograms of the mass, radius, Umax and 
Ai/ of the simulated CHeB stars with a = 2.4. An increase in a re- 
sults in a movement of the peak location of the radius distribution 
towards a lower value, and leads to a movement towards a higher 
value for the peak locations of the Vmax and Av distributions. Fur- 
thermore, Fig.|9]represents the distributions of the fmax and Av as 
a function of stellar age. Increasing the mixing-length parameter a 
leads to the upper boundary of the distributions moving up slightly. 



4 DISCUSSIONS AND CONCLUSIONS 

Several points should be kept in mind when comparing the theo- 
retical distributions of Vmax and Au with the observations. Firstly, 
the scaling equation (|TJ is obtained from the ideas that i^max is pro- 

— 1/2 

i iortional to the aco ustic cutoff frequency Vac and i^ac oa gT^ , / 
Brown et al.ll9 91 ). where g is the stellar surface gravity and Teff 
is the effective temperature. The Umax given by e quation l[T) is not 
necessarily accurate JGillilandl2008l ; IStello et al.ll2009al) . For ex- 
ample, the ratio of Uac of stellar models to that calculated from the 
global p arameters of models is always below unity ( see the Fig- 
ure 8 of IStello et al.l 12009a|)). However, the scaling equation {2} 
is very a ccurate. The uncert ainty of equation ([2} is within a few 
per cent dStello et al.ll2009al la). In addition, scaling equations (TJ 
and ^ agree within a few per cent with stellar model calculations 
for cool models (T^ff < 6400 K) JStello et alj[2009ah . A few per 
cent deviations in equations ([TJ and ^ cannot obviously change 
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Figure 9. The distributions of the Umax and Ai^ of the CHeB stars of sim- 
ulated SSP with a = 2.4 as a function of stellar age. 



the peak locations of Umax and Au in our simulations. Hence, the 
deviations cannot affect our basic results. 

The giant stars observed by the CoRoT may contain the first 
giant branch (FGB) stars. Our simulations show that the distribu- 
tions of Umax and Au of the FGB stars are uniform. Therefore, 
the Umax and Ai^ of the FGB stars do not af fect the peak locations. 

Comparing with the lMiglio et alj ( l2009h 's simulation, we used 
the same IMF and age-metallicity relation. In our simulation, the 
SFR is a constan t over the pas t 15 Gy r. But it is a constant over the 
past 9 Gyr in the lMiglio et al.l ( l2009n 's simulation. This difference 
could not affect the simulated results because the oldest population 
are impeded into the CHeB stage by mass loss. For example, for ry 
= 1.0, the simulated population do not contain the stars older than 
about 10 Gyr. The f requencies of our sim ulated peak locations are 
lower than those of iMiglio et alj ( 120090 . This may be caused by 
mixing-length parameter. 

For different rj, the CHeB stars with age < 2 Gyr lose only a 
negligible amount of mass from their surface before they evolved 
into CHeB stage. Consequently, the effect of the mass loss on the 
evolution of these stars is insignificant. The distributions of Umax 
and Au of these stars are almost not affected by Reimers mass 
loss. However, for old population, the bigger the rj, the more the 
mass loss. For a big rj, stars can lose an appreciable, but from star 
to star different, amount of mass from their surface before the he- 
lium flash. Then the stars move to the left of the HR diagram with 
slightly decreasing luminosity after the flash, and the mean density 
of the stars is enhanced. The increase in the mean density leads 
to an increase in the Umax and Au of these stars. However, the 
change in the mean density is different from star to star. The mass 
of the stars evolved into CHeB at the same age is approximate, but 
the extent of central helium burning is different. The stars around 
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the ZAHB have a smaller radius than those approaching the AGB. 
The fractional change caused by mass loss in the radius of the stars 
around the ZAHB is larger, and the change in the mean density of 
these stars is also larger than that of those approaching the AGB. 
The stars which are close to the ZAHB are located around the up- 
per boundary of the distributions of v max and Av with age, while 
those approaching the AGB are located around the lower boundary. 
Therefore, the changes of the upper and lower boundary shown in 
Fig. |5] are different. However, for the middle age population, the 
stars lose only a little amount of mass from their surface before the 
helium flash, which hardly affect the radius and luminosity of the 
stars after the flash. Because the mass decrease slightly but the ra- 
dius is almost not changed, the mass loss results in that the Umax 
and /S.v of these stars decrease slightly. 

The Reimers mass loss mainly affects the old population. The 
mass loss leads to an increase in the Vmax and Av of old popu- 
lation. However, a high mass-loss rate can impede the low-mass 
low-metal stars evolving into CHeB stage. If the Reimers mass loss 
is very efficient during the red-giant branch, the Vmax and Ajy of 
CHeB stars would not be observed in old clusters. The asteroseis- 
mical observation on old clusters may provide a help to constrain 
the mass-loss rate. 

For 77 = 0.5, the value of Vmax and Ai/ of population with 
age > 2 Gyr is almost located in the range of 10-30 ^Hz and 1-4 
yuHz, and gathers about 25 /iHz and 3 /xHz, respectively. Therefore 
increasing or decreasing the number of stars at a certain age > 2 
Gyr cannot affect the peak locations of the Vmax and Av. So the 
peak locations are not sensitive to the SFR and not sensitive to 
whether population contains old stars. For r) = 1.0, increasing the 
SFR of stars with age between 7-9.5 Gyr can increase the stars with 
I'max > 25 /iHz and Av > 3 /iHz. But even enhancing the SFR to 
several times, the peak locations are not affected. 

In BSP, some CHeB stars can lose a little mass, but some stars 
can accrete a little mass by the weak binary interactions. The effect 
of a slight mass change before the helium flash on the luminosity 
and radius of CHeB stars after the flash is negligible. Therefore 
the Umax and Av of the CHeB stars which lost a little bit of mass 
decrease; whereas those of the CHeB stars which accreted a little 
mass increase. Hence the distributions of Vmax and Ai/ of CHeB 
stars cannot be affected. However, for the strong binary interac- 
tions, on the one hand, some stars which have lost an appreciable 
amount of mass from their surface would move to the left of the 
H-R diagram with slightly decreasing luminosity, at the same time, 
the mean density of these stars would increase. Consequently the 
Vmax and Av of these stars increase. On the other hand, some stars 
with age > 2 Gyr which have accreted a considerable amount of 
mass would become like the stars with age < 2 Gyr, hence their 
I'max and Au can increase obviously too. However, the fraction of 
these interactive binary stars appearing in our simulated CHeB stars 
is very small. Therefore, although the binary interactions such as 
mass transfer and mass accretion can affect the Vmax and Au of 
CHeB stars undergoing a mass accretion or mass loss, the effect 
of binary interactions on the distributions of the Vmax and Av of 
CHeB stars is not significant. 

An increase in the mixing-length parameter a mainly leads to 
a decrease in the radius of all CHeB stars after the helium flash. 
Thus the i'max and Av of CHeB stars increase with a. Conse- 
quently, the peak locations of the Vmax and Ai/ can be affected 
by the a. For the middle-age population at the age of about 2 to 5 
Gyr, the peak locations are more sensitive to the a than mass-loss 
rate, SFR and metal abundance. Thus the asteroseismical observa- 



tion on the middle age clusters may provide a help to constrain the 
mixing-length parameter. 

The mass of simulated CHeB stars is mainly located in 1-2 
Mq. For the CHeB stars with the mixing-length parameter a = 
2.0, the most of them have a radius in the range of 11-14 Rq. 
Therefore, even the mass distribution is uniform, the stars have an 
approximate mean density. Hence they have an approximate i^max 
and Au. Consequently, there is a dominant peak in the distributions 
of the Vmax and Av. 

A high Reimers mass loss can lead to an increase in the Vmax 
and Av of old population, and impedes the low-mass low-metal 
stars evolving into CHeB stage. On the contrary, it results in a 
very small decrease in those of middle-age population. However, 
the mass loss scarcely affect the Vmax and Av of young popula- 
tion. The effect of the Reimers mass loss on the peak locations of 
Vmax and Av of the CHeB stars is not significant unless the mass- 
loss rate is very high. The effect of star formation rate and binary 
interactions on the peak locations of Vmax and Av of the CHeB 
stars is also not significant. The dominant peak of Vmax and Av 
is due to the fact that most of CHeB stars have an approximate ra- 
dius. The radius can be affected by the mixing-length parameter. 
The peak location also can, thus, be affected by the mixing-length 
parameter. 
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